function c = makesynthwaves(n)
   
   % Add N traces of synthetic data to a correlation object. Alters properties: waves, trig,
   % start, Fs.
   
   % Author: Michael West, Geophysical Institute, Univ. of Alaska Fairbanks
   % $Date$
   % $Revision$
   
   
   % MAKE SINGLE SIGNAL
   t = [-5:.01:14.99]';
   env_orig = exp(-1*(t-5).^2/ (2*3^2) );
   sprep = 0.2*sin(t+1) + .3*sin((t-3)/0.3) + 0.4*cos(t/0.8) + 0.2*cos((t-1.7)/0.2);
   %sprep = sin(t);
   %sprep = 0.5*sin(t) + 0.5*sin(t/.5);
   s = sprep .* env_orig;
   
   
   % MAKE SUITE OF SIMILAR EVENT
   %n = 150; %no. of events
   w = s * ones(1,n) + .01*(rand(length(s),n)-.5); % add time offset
   aa = .6;    % adjust to add randomness
   for i = 1:n
      w(:,i) = w(:,i) + aa*rand(1)*sin((t-rand(1))/rand(1));
   end;
   
   
   % VARY START TIMES
   for i = 1:n
      bump = round(300*(rand(1)));
      w(:,i) = w([ length(t)-bump:length(t) 1:(length(t)-1)-bump],i);
   end;
   
   
   % VARY AMPLITUDES
   for i = 1:n
      w(:,i) = rand * 100 * w(:,i);
   end;
   
   
   % DEMEAN
   w = w - ones(size(w,1),1)*mean(w);
   
   
   % MAKE v0 CORRELATION
   d.trig = 732604 + rand(n,1);
   d.start = d.trig - (5 + 0.2*rand(n,1))/86400;
   d.Fs = 100;
   d.w = w;
   
   
   % CREATE v1 CORRELATION OBJECT
   c.W = waveform;
   for i = 1:length(d.trig)
      w = waveform;
      w = set(w,'Station','UNKN');
      w = set(w,'Channel','EHZ');
      w = set(w,'Start',d.start(i));
      w = set(w,'Fs',d.Fs);
      w = set(w,'Data',d.w(:,i));
      c.W(i) = w;
   end;
   c.W = reshape(c.W,length(c.W),1);
   c.trig = d.trig;
   
end
